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The surface partition of large clusters is studied analytically within a frame- work of the "Hills and 
Dales Model". Three formulations are solved exactly by using the Laplace- Fourier transformation 
method. In the limit of small amplitude deformations, the "Hills and Dales Model" gives upper 
and lower bounds for the surface entropy coefficient of large clusters. A comparison with the 2- and 
3-dimensional Ising model surface entropy coefficients is made. 



I. INTRODUCTION 

The surface entropy of large clusters was introduced 
by Fisher in his droplet model (FDM) Ij. From the 
study of the combinatorics of clusters, Fisher postulated 
that the leading contribution to the surface entropy is 
proportional to the surface S, i.e. ujS (in dimensionless 
units). The coefficient w is a ratio of the surface energy 
coefficient cro(Tc) per one constituent taken at the crit- 
ical temperature Tc. The surface entropy was studied 
recently in our paper There we developed the "Hills 
and Dales Model" (HDM) which is a statistical model 
of surface deformations that obeys the volume conserva- 
tion of considered clusters. Using the novel mathematical 
method, the Laplace-Fourier transform Q, we were able 
to find the grand canonical surface partition (GCSP) of 
the HDM analytically. For vanishing deformations we ob- 
tained the upper limit for the surface entropy coefficient 
u! of large clusters to be w ~ 1.06009 (in dimensionless 
units), i.e. about 6 % larger than Fisher's postulate. 

In the grand canonical formulation the cluster vol- 
ume is conserved in average, but to apply the HDM to 
small and finite clusters it is necessary to consider an ex- 
act volume conservation using the canonical formulation. 
Therefore, in present paper we consider the canonical for- 
malism for the HDM and obtain the infimum for the sur- 
face entropy of finite and large clusters. For the limit of 
vanishing deformations we also introduce the semi-grand 
canonical ensemble which occupies an intermediate place 
between the grand canonical and canonical surface en- 
sembles. With the help of the Laplace-Fourier transform 
technique Q the canonical surface partition (CSP) and 
the semi-canonical surface partition (SGCSP) are evalu- 
ated exactly for any volume of cluster. For large clusters 
the leading contribution and its corrections are found an- 
alytically for the CSP and SGCSP. The obtained values 
for the cj-coefficient are compared with the corresponding 
values for 2- and 3-dimensional Ising models for different 
lattice geometries. It is shown that u! values of all 2- 
and 3-dimensional Ising models lie between the supre- 
mum and infimum of the HDM. 

The paper is organized as follows. In Sect. II we for- 
mulate three ensembles for surface deformations within 
the HDM frame work and solve them analytically by the 
Laplace-Fourier transform technique. Sect. HI is devoted 



to the analysis of the isochoric ensemble singularities and 
to the derivation of the upper estimates of the surface 
entropy coefficient. The lower estimates for the surface 
entropy coefficients are found and compared to the cor- 
responding 2- and 3-dimensional Ising lattice values in 
Sect. IV. The conclusions are formulated in Sect. V. 



II. HILLS AND DALES MODEL 

The HDM is a statistical model of surface deforma- 
tions. We impose a necessary constraint that the defor- 
mations should conserve the total volume of the cluster 
of A-constituents. As in our previous paper , the main 
interest is focused on the deformations of vanishing am- 
plitudes. This is sufficient to find both an absolute supre- 
mum and absolute infimum for the w-coefficient of the 
HDM. In this case the shape of the deformation cannot 
be important to our result, so we can choose the regular 
one. For this reason we shall consider cylindrical defor- 
mations of positive height > (hills) and negative 
height —hk (dales), with fc-constituents at the base. For 
simplicity it is assumed that the top (bottom) of the hill 
(dale) has the same shape as the surface of the original 
cluster of A-constituents. We also assume that: (i) the 
statistical weight of deformations exp {—ao\ASk\/ si/T) 
is given by the Boltzmann factor due to the change of 
the surface lASfel in units of the surface per constituent 
si; (ii) all hills of heights hk < Hk {Hk is the maximal 
height of the hill with fc-constituents at the base) have 
the same probability dhk / Hk besides the statistical one; 
(iii) assumptions (i) and (ii) are valid for the dales. 

These assumptions are not too restrictive and allow us 
to simplify the analysis and find the one-particle statisti- 
cal partition of the deformation of the fc-constituent base 
as a convolution of two probabilities discussed above: 
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(1) 



where upper (lower) sign corresponds to hills (dales). 
Here Pk is the perimeter of the cylinder base. Our next 
step is to find a geometrical partition (degeneracy factor) 
or the number of ways to place the center of a given defor- 
mation on the surface of the A-constituent cluster which 
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is occupied by the set of {nf^ — 0, 1, 2, ...} deformations 
of the Z-constituents base. 

For the grand canonical surface partition (GCSP) the 
desired geometrical partition can be given in the excluded 
volume approximation j^j as 



Gnc = 



K„ 



Sa- E k{n+ + n^^ ) si 
fe=i 
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(2) 



where sik is the area occupied by the deformation of k- 
constituent base (fc = 1,2,...), Sa is the full surface of 
the cluster, and Kj^axi^A) is the A-dependent size of the 
maximal allowed base on the cluster. It is clearly seen 
now that the first multiplier in the right hand side (r.h.s.) 
of (0) corresponds to the available surface to place the 
center of each of {n^} deformations that exist on the 
cluster surface. It is necessary to impose the condition 
Qgc > which ensures that the deformations do not over- 
lap. Equation (0) is the van der Waals excluded volume 
approximation usually used in statistical mechanics at 
low particle densities 0, B IS ^ and can be derived 
for objects of different sizes in the spirit of Ref. [l^. 

According to Eq. iQJ the statistical partition for the 
hill with a fc-constituent base matches that of the dale, 
■ - -+ ~ and, therefore, the GCSP 
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corresponds to the conserved (on average) volume of the 
cluster because the probabilities of hill and dale of the 
same base are identical. The 6(si5gc)-function in © 
ensures that only configurations with positive value of 
the free surface of cluster are taken into account, but 
this makes the calculation of the GCSP very difficult. 

For small and finite clusters we have to impose a 
more strict constraint of the exact volume conservation 
of cluster. This can be done within the canonical en- 
semble assuming that the number of the hills n'^ of the 
/c-constituent base is always identical to the number of 



corresponding dales, i.e. 



rifc. Then the 



canonical geometrical partition can be cast as follows 



Gc = 
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(4) 



where the factor two in the denominator of the right hand 
side (r.h.s.) of Q accounts for the fact that it is necessary 
to place simultaneously the centers of two fc-constituent 
base deformations (hill and dale) out of 2nk on the sur- 
face of cluster. Using the geometrical partition |@J, one 
can write the partition function of canonical ensemble, 
the CSP, as follows 
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As in the case of the GCSP, the 0(2si^c)-function in the 
CSP ensures that only configurations with positive value 
of the free surface of cluster are accounted for, but this 
constraint makes the calculation of the partition (jSJ very 
difficult. 

An additional problem in evaluating the partitions 
and jSJl appears due to the explicit Sa dependence of the 
maximal base of deformations via Kmax{SA) because in 
this case the standard method to deal with the excluded 

volume partitions, the usual Laplace transform ^ 

[Tol | in Sa, cannot be applied. However, as shown in 
the GCSP © can be solved analytically with the help 
of the Laplace- Fourier technique 0. The latter employs 
the identity 
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which is based on the Fourier representation of the Dirac 
(5-function. Similar to the GCSP, the representation © 
allows us to decouple the additional S'^-dependencc in 
Kmax{SA) of the CSP and reduce it to the exponential 
one, which can be integrated by the Laplace transform 
ill 
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After changing the integration variable Sa 

Sa — 2 E kukSi, the constraint of 0-function has 
fc=i 

disappeared. Next all Uk were summed independently 
leading to the exponential function. Now the integration 
over S' in (0 can be done giving the canonical isochoric 
partition 
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where the function J^c(^, A) is defined as follows 

= (f) ..+ 
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The representation ((SJ is generic - it is also valid for the 
GCSP, if the canonical function @ is replaced by the 
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grand canonical one 
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Before making the inverse Laplace transform and 
studying the structure of singularities of the functions 
© and pOfl . it is necessary to discuss one more ensem- 
ble for the surface deformations which hereafter will be 
called as the semi-grand canonical surface partition. 

This ensemble occupies an intermediate position be- 
tween the canonical and grand canonical formulations. It 
corresponds to the case, when the hills and dales of the 
same base are considered to be indistinguishable. For 
that we would like to explore the fact that according to 
the statistical probabilities of hills and dales of the 
same base are equal. Then for the infinitesimally small 
amplitudes of deformations the volume conservation con- 
straint is fulfilled trivially. For finite amplitudes of defor- 
mations this ensemble corresponds to the case of clusters 
of intermediate size, when its surface is sufficiently filled 
with deformations that it is impossible to add one more 
deformation of large base in order to conserve the cluster 
volume, but it is possible to place many deformations of 
a much smaller (or even smallest) base. 

For this case the geometrical factor reads as 



k=l 
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and the SGCSP has the following form 
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It is easy to show that using the Laplace-Fourier trans- 
form technique 3] the SGCSP (|12|l can be transformed 
into the generic representation (jS)) for the function 



k=i *i 



k siX 
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By construction the equations H12I) and 11311 arc less fun- 
damental than the corresponding grand canonical and 
canonical functions. 



III. ANALYSIS OF SINGULARITIES 

To study the structure of singularities of the isochoric 
partition JHJl, it is necessary to make the inverse Laplace 
transform {a G {gc, sg, c}): 
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where the contour integral in A is reduced to the sum over 
the residues of all singular points \ = \n + irj with n = 
0, 1, 2, .., since this contour in the complex A-plane obeys 
the inequality x > max(i?e{A„}). Now all integrations 
in (|14|l can be done, and all three surface partitions {a G 
{gc, sg, c}) can be written as 
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i.e. the double integral in H14() simply reduces to the 
substitution ^ — > Sa in the sum over singularities. This 
remarkable answer for all three surface partitions is a 
partial example of the general theorem on the Laplace- 
Fourier transformation properties proved in 

The simple poles in (|14|l are defined by the condition 
An = Ta{SA, A„) and for each ensemble the latter can be 
cast as a system of two coupled transcendental equations 

E C e-^'''"Cos(/;^fc), (16) 

k=l 
K^uASa) 

i: = - E '/'fe e-'=«"sin(/:fc), (17) 
fe=i 

for dimcnsionless variables defined as i?" = sii?e(A„) 
and = si/m(A„) for the GCSP and SGCSP, and as 
= 2sii?e(A„) and = 2si/m(An) for the CSP. Here 
the function t/)^ is given by the expression 



'-k ' 



+ Zj, , for a ~ gc, 
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To this point Eqs. and (dJ) are general 

and can be used for particular models which specify the 
height of hills and depth of dales. But it is possible to 
give both the upper and lower estimates for all three 
partition functions of large clusters, and even to esti- 
mate corrections for finite and small clusters. For the 
upper estimate let us consider the real root (i?g ; Iq = 0) 
of these equations. It is sufficient to consider the limit 
Kmax{SA) OO, bccausc for /" = /q = the r.h.s. of 
(|16|l is a monotonously increasing function of KmaxiSA)- 



Since 



are the monotonously decreasing func- 



tions of Hk, the maximal value of the r.h.s. of (|16|) cor- 
responds to the limit of infinitesimally small amplitudes 
of deformations, Hk ^ ^ zt = zZ = 1 . Under these 
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conditions Eq. p7|) for I" = 1^ = becomes an identity 
and Eq. (|16|l acquires the form 



(19) 



k=l 



where BS" = 2 and B" ^ B'9 ^ 1. Therefore, the real 
roots of and the corresponding surface entropy co- 



efficients ojj} are as fohows 



Hq - 
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1.060090, 

0.806466 , 
0.806466 , 



for a = gc, 
for a — c, (20) 
for a = sg . 



The resuhs (|20|) correspond to the upper estimate for 
the surface partitions because for /" ^ defined by lfT7|) 
the inequafity cos{I"k) < 1 cannot become the equahty 
for all values of k simultaneously. Then it follows that 
the real root of ((TB|l obeys the inequality Rq > i?">o- 
The last result means that in the limit of infinite cluster, 
S'a — * oo, all surface partitions H15|) are represented by 
the farthest right singularity among all simple poles {A„} 



max 
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1 



J^g(i?.g+i3°) 



= 9ae 



, ,a Sa 
OJjj — 



(21) 



where the geometrical degeneracy prefactor g^ is defined 
as follows: gg^ « 0.38139 and g^ = g^g « 0.407025. Thus, 
the geometrical factor of the leading term for all three 
models is practically the same. 

Remarkably the result H21|l is model independent. This 
is a consequence of the vanishing deformations limit in 
which all model specific parameters vanish. The second 
remarkable fact is that Eq. (|21|l is valid for any self-non- 
intersecting surfaces of cluster. This is so because both 
the shape and dimensionality of the cluster under con- 
sideration do not enter in our equations explicitly. For 
our analysis of the HDM surface partitions it was suffi- 
cient to require that the cluster surface together with de- 
formations is a regular surface without self-intersections. 
Therefore, for vanishing deformations the latter means 
that Eq. 1)2 l|l should be valid for any self-non-intcrsecting 
surfaces. 

For large, but finite clusters it is necessary to take into 
account not only the farthest right singularity Aq in (|15|l . 
but all other roots with positive real part i?">o > 0. The 
analysis presented in Appendix A shows that besides the 
opposite signs there are two branches of solutions, 
and " , for the same n > 1 value: 
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FIG. 1: The first quadrant of the complex plane siA„ = 
Rn +iln shows the roots of the system of Eqs. (Ilfcil and (1171 . 
The symbols represent the two branches and of the 
roots for the upper estimate of three surface partitions. The 
curve is defined by the approximation given by 1221 and 12311 
(see text for more details). 



curve parametrized by functions /" and i?" taken from 
Eqs. and (|23J), respectively. From Eq. and 

Fig. 1 it is clear that for the GCSP the largest real 
part w 0.0582 is about 18 times smaller than R^'^, 
whereas for the CSP and SGCSP the real part i?f = i?f 
of the first most right complex root of Eqs. H16|) and H17|) 
is about 63.6 times smaller than i?§ = Rq^ . Therefore, 
for a cluster of a few constituents the correction to the 
leading term 121|l is exponentially small for all considered 
partitions. Using the approximations (|22|l and (|23|l . for 
n > 2 one can estimate the upper limit of the {R!^;I^ 
root contribution into Eq. I|15|l 



Sa 



< e s^'^n-s, I . (24) 



This result shows that for all three considered partition 
the total contribution of all complex poles in (|15|l is neg- 
ligibly small compared to the leading term H21(l for a 
cluster of a few constituents or more. 



IV. THE SURFACE ENTROPY COEFFICIENTS 



RZ 



(B^ 



(23) 



The exact solutions (i?";/" ) for n > 1 which have the 
largest real part are shown in Fig. 1 together with the 



To complete our analysis of the limit of vanishing de- 
formations we would like to find the lower estimate for 
the GCSP, CSP and SGCSP for large clusters. This es- 
timate corresponds to the absence of all other deforma- 
tions except for those of smallest base. In other words. 
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one has to substitute K^axiSA) = 1 in all corresponding 
expressions. Then equations 1)16(1 and l|17|) . respectively 
become 



e-«" cos(7«) , 
-0? sin(/;^) . 



(25) 
(26) 



Similar to the previous consideration, the leading term of 
the lower estimate for the surface partitions ((15(1 is given 
by the real root {R^;I^ = 0) of the system ^ and it^ 
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2luI = 2min{w^} 
Lo'^jf = minjaj'^f } 
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for a 
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c, (27) 
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FIG. 2: The second quadrant of the complex plane siA,i = 
Rn + iln shows the complex roots of the system of Eqs. 1(25(1 
and 1261 1 with the largest real parts. The circles and squares 
represent the roots for the lower estimate of the GCSP and 
CSP(SGCSP), respectively 



Partition 


maxjcj"} 


min{cj"} 


GCSP 


1.060090 


0.852606 


SGCSP 


0.806466 


0.567143 


CSP 


0.403233 


0.283572 



Table I. The maximal and minimal values of the oj-coefiicient 
for three statistical partitions of the HDM. 

The ^-coefficients for upper and lower estimates of all 
three surface partitions are summarized in the Table I. 
A comparison with the corresponding coefficient for liq- 
uids should be made with care because various contribu- 
tions to the surface tension, i.e., eigen surface tension of 
the liquid drop, the geometrical degeneracy factor (sur- 
face partition), and the part induced by interaction be- 
tween clusters, are not exactly known. Therefore, even 
the linear temperature dependence of the surface tension 
cr(r) = ao{Tc - T)/Tc due to Fisher applied to a nu- 
clear fiquid ( (To « 18 MeV; « 18 MeV Q) may be used 
to estimate the w-coefficient, if both the eigen surface ten- 
sion and the interaction induced one are non-increasing 
functions of temperature. Under these assumptions one 
can get the following inequality for nuclear liquid 



UJ„ 



1.060090, 



(28) 



i.e. the upper estimate for the GCSP, indeed, provides 
the upper limit for surface partition of nuclear matter. 

A similar analysis for real liquids is difhcult because 
of complicated temperature dependence of the surface 
tension. Therefore, we would like to compare the uj- 
coefficients from Table I with the w-coefficients for the 
large spin clusters of various 2- and 3-dimensional Ising 
models, which are listed in the Tables II and III, respec- 
tively [13 ■ Such a comparison can be made because the 
surface entropy of large spin clusters on the Ising lattices 
are similar to the considered surface partitions 1(15(1 [T^ . 

The (jj-coefhcient for the d-dimensional Ising model is 
defined as the energy 2J required to flip a given spin 
interacting with its g-neighbors to opposite direction per 
(d— l)-dimensional surface divided by the value of critical 
temperature 



1± 



(29) 



Again, as in case of upper estimates one can show that 
the real root (i?g ;Iq =0) approximates well the lower 
estimate for the partition function for a system of a few 
constituents. In fact, each of three surface partitions has 
only a single root with positive real part which coincides 
with (i?Q ; Iq ~ 0). In Fig. 2 a few complex roots of Eqs. 
((25(1 and l(26() with the largest real parts are shown. Since 
all these roots have negative real part, they generate an 
exponentially small contribution to the lower estimate of 
surface partition for a system of a few constituents. 



Here q is the coordination number for the lattice, and 
J denotes the coupling constant of the model. A com- 
parison of the Tables I - III shows that all lattice ujLat- 
coefficients, indeed, lie between the upper estimates for 
canonical and grand canonical surface partitions 



0.403233 



< UJLat < W 



U 



1.060090, (30) 



i.e. uJij and Luf^ are the infimum and supremum for 2- 
and 3-dimensional Ising models, respectively. 
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Lattice type 


^Lat — ^ 


Honeycomb 


0.987718 


Kagome 


0.933132 


Square 


0.881374 


Triangular 


0.823960 


Diamond 


0.739640 



Table II. The values of the ojLat-coefScient for various 2- 
dimensional Ising models. For more details see the text. 



Lattice type 


(^Lat — 


Simple cubic 


0.44342 


Body-centered cubic 


0.41989 


Face-centered cubic 


0.40840 



Table III. The values of the tJLat-coefEcient for various 3- 
dimensional Ising models. 

The HDM partitions do not have an explicit depen- 
dence on the dimension of the surface, but a compar- 
ison of the HDM and Ising model w-coefRcients shows 
that the HDM ensembles seem to posses some sort of 
internal dimension: the GCSP is close to honeycomb, 
kagome or square lattices, whereas the SGCSP is simi- 
lar to triangular and diamond lattices, and the maxjo;} 
of the CSP is closer to the 3-dimensional Ising models. 
In some cases the agreement with the lattice data is re- 



markable 



coincides with the arithmetical average 



of the w-coefHcients for square and triangular lattices up 
to a fifth digit, but in most cases the values agree within 
a few per cent. The latter is not surprising because the 
HDM estimates the surface entropy of a single cluster, 
whereas on the lattice the spin clusters do interact with 
each other and this, of course, changes the surface tension 
and, consequently, affects the value of critical tempera- 
ture. It is remarkable that so oversimplified estimates of 
the surface partitions for a single large cluster reasonably 
approximate the w-coefficients for 2- and 3-dimensional 
Ising models. 

It would be interesting to check whether the lower es- 
timate of the CSP w£ « 0.283572 is an infimum for the 
Ising lattices of higher dimensions d > 3. If this is the 
case, then we can give an upper limit for the critical tem- 
perature of those lattices using Eq. H29|l 



J 



< 



3.5264 , 
d 



(31) 



On the other hand, the lower estimate for the critical 
temperature of Ising lattices, ^ > ^/e ^ , is provided by 
the supremum of the w-coefficients of surface partitions. 



V. CONCLUSIONS 

We have formulated the grand canonical and canoni- 
cal partitions of surface deformations in the frame work 
of the HDM. Both partitions conserve the volume of the 
deformed cluster and take into account all surface de- 
formations with non-negative value of the free surface of 
this cluster. The grand canonical surface partition con- 
serves the cluster volume in average, whereas in canonical 
formulation it is conserved exactly. These partitions are 
solved exactly for an arbitrary (finite or infinite) size of 
largest deformation by the Laplace-Fourier transforma- 
tion technique and the the general analytical expression 
(|15|l for these partitions in terms of the set of isochoric 
ensemble singularities are derived. 

Similarly we introduced and solved a special ensem- 
ble, a semi-grand canonical partition, which obeys all 
constraints discussed above in the limit of vanishing de- 
formations and occupies an intermediate place between 
the grand canonical and canonical ensembles. 

Then we considered the limit of vanishing deformations 
for all three surface partitions, and obtained the upper 
and lower estimates for the surface entropy for each of 
these partitions. A comparison of obtained w-coefiicients 
for surface partitions with the corresponding coefficients 
for the large spin clusters of 2- and 3-dimensional Ising 
models shows that the upper estimate of the GCSP is a 
supremum, whereas the upper estimate of the CSP is an 
infimum for considered lattices. The question of the Ising 
lattice w-coefficients for higher dimensions is discussed. 

The developed formalism is rather general and, there- 
fore, may be applied to the surface deformations of any 
kind of clusters, if the underlying mechanism of the sur- 
face deformations is given. 

Acknowledgments. This work was supported by the 
US Department of Energy. 



VI. APPENDIX A 

For large, but finite clusters it is necessary to take into 
account not only the farthest right singularity Aq in l|15|l. 
but all other roots of Eqs. (|16|l and H17|) which have 
positive real part i?">o > 0- In this case for each i?">o 
there are two roots ±7" of (|17|l because the partition 
function (|15|) is real by definition. The roots of Eqs. 
(|16|l and (|17|) with largest real part are insensitive to 
the large values of K^axiSA), therefore, it is sufficient 
to keep Kmax{SA) — > oo. Then for limit of vanishing 
amplitude of deformations Eqs. (|16|l and (|17|l can be, 
respectively, rewritten as 



e^" cos(/,';')-l 



= - e^<^ sin(4") 



(32) 
(33) 
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After some algebra the system of (|32|l and \3'6\i can be 
reduced to a single equation for R" 



_B°(B°+2_R°) 



a\2 



1/2 



cosh R" 



2 i?" 



sinhi?^, (34) 



two branches of solutions, J""*" and for the same 

n > 1 value. Expanding both sides of for i?" ^ 1 
and keeping the leading terms, one obtains H22|) and ^'2.'A\ . 
In Fig. 1 this approximation is compared with a few 
exact solutions /^) for n > 1 which have the largest 
real part. 



and the quadrature = y^ ^"<f^°„+_^^^"^ - [KY- The 
analysis shows that besides the opposite signs there are 
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